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The formulas of field strength over an inhomogeneous spherical earth are obtained on 
the conditions that (I) the radius and the electrical properties of the earth's surface dis- 
continuously change several times along the wave path, or (II) the surface of terrain arbi- 
trarily changes in height along the wave path, but it is still smooth everywhere and the 
radius of curvature is sufficiently large as compared with the wavelength. 

The case (I) is considered to be more general than those of mixed pal lis on a smool h 
earth, because the latter can be seen as special cases of the former. The case (II) cor- 
responds to the case of multiple diffraction of radio waves by several mountains having 
finite radii of curvature. In Doth cases, the unified formulas of field strength are obtained 
in the form of a multiple residue series, which is reduced to the ordinary Van der Pol and 
Bremmer formula in the special case of homogeneous ground. 

The convergence of series of the formulas is very good when the propagation distance 
on every section of the inhomogeneous ground is long enough or the diffraction loss is large 
enough, and is poor when any one of these distances is so short that the section is effectively 
seen as a flat plane, or the diffraction loss on the section is very small. In these cases, the 
flat earth or other approximations can be used, and several supplementary formulas are 
prepared for cases of poor convergence. Several special applications are given. 

Introduction 

The problems of mixed paths over a smooth earth have been investigated by many authors, 
especially in the case of a Sal earth, and many equations and approximations have been estab- 
lished according to the given situation. However, the mixed paths in these eases mean tlie 
radio wave propagation path over a smooth ground, in which several sections having different 
electrical properties are included, and thus changes in height arc not included. 

In part I of this paper, these changes in height are taken into account on the assumption 
that the ground surface changes discontinuously along the wave path, as illustrated in figure 3. 
A series of papers has been established on the theory of propagation over terrain of this model 
[Furutsu, 1957a, 1957b, 1959a], and the unified formula of field strength was obtained in the form 
of a multiple residue series, which reduces to that of mixed paths in the special case of smooth 
earth. This model of terrain form may not be suitable for the range of very high frequency 
where the radius of curvature of the terrain becomes sufficiently large as compared with the 
wavelength, and it also may not be suitable for a completely irregular terrain where some 
statistical approach would be more appropriate. In the former ease, the problem could be 
treated as a multiple diffraction by hills or mountains having large radii of curvature, in which 
the earth's surface could be included, as illustrated in figure 23. Also in this ease, the field 
strength can be obtained in a unified form of multiple residue series [Furutsu, 1956], and part II 
of this paper is devoted to this subject. 

The purposes of this paper are the survey of these formulas and the several possible appli- 
cations to the practically important cases, such as the evaluations of the effects on the ground 
radio wave of a ridge, a cliff, and a bluff at a close line (part I) or those of several hills or moun- 
tains, taking into account their curvatures if necessary (part II). 



1 This work was sponsored by the Air Force Cambridge Research Laboratories, Office of Aerospace Research (USAF), Bedford, Massachusetts. 
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Part I. Radio Wave Propagation Over Inhomogeneous Spherical Earth 

1. Equation Formulation 

Solving- Maxwell's equation in terms of the vector potential, A, the equation in the Car- 
tesian coordinate system (x) = (x, y, z) becomes (the earth surface is assumed to be flat for the 
time being) 

l(yk- 2 v)+l]A(x) = -k-^J(x), F=coWo(e + <^), v= (^ 7 1,'^} (L1) 

Here, the time factor e im is omitted and /is the current density of the external system. Other 
notations are those ordinarily used. 

In the case of the vertical dipole excitation at the arbitrary point x', we may first presuppose 
that the horizontally polarized wave will not be induced by the terrain of the form assumed 
and thus, using the 2-axis in the vertical direction from the ground, put 

A(x) = (0, 0, *(*,«')), k- 2 fjLoI(x) = (0, 0, B(x-x')). (1.2) 

This presupposition will be justified later. Thus, the solution can be represented by the Green 
function defined by the equation 

[(vk- 2 v) + l}Hx,x') = -5(x-x'). (1.3) 

Here, the boundary condition on every horizontal surface made by the medium discontinuities 
is found to be the continuity condition of 

m, (*- f J» a-*) 

where d/dn=d/dz and the earth's surface is flat in this case. 

Even when the curvature of the earth's surface is taken into account, eq (1.3) is still valid 
in a small but practically sufficient domain of the earth's surface; it is sufficient to assume, for 
any small line element having the components (dx, dy, dz), the length ds given by the metrics 

ds 2 =^z/a) 2 (dx 2 +dy 2 )+dz 2 . 

Here, the coordinate system is taken so that the earth's surface is given by the surface z=a, a 
being the earth's radius. Further, the boundary condition (1.4) is also valid on the condition 
that the earth's radius is sufficiently large as compared with the wavelength. 

In spite of the fact that eqs (1.3) and (1.4) are valid only in the case where the medium 
changes with the ^-coordinate but not with the x, ^/-coordinates, we are going to treat the case 
where the height and electrical properties of the earth's surface take spatially discontinuous 
values, as is illustrated in figure 3. Indeed, from the definition of \p, it is readily proven that 
the boundary condition (1.4) does not hold on the vertical boundary surfaces of medium dis- 
continuity. However, as is seen in the following, this fact does not give any serious effect 
to the result of the assumption that eq (1.4) also holds on the vertical boundary surfaces. 

For any continuous functions of \p' and ^", the Green theorem can be given in the form 

(%. i) 

f ^ / Figure 1. The integration domain for eq {1.5). 
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Here, the left side is the volume integral in the space 2 and 

w(s),i"{s)]=-w{s),ris)] 



(1.5) 



(1.6) 



Here, d/d/^ is the inward normal differentiation on the surface s of the arbitrary space 2. 

It is especially to be noted that the lemma (1.5; holds even when k takes discontinuous 
values across the surfaces contained in the space 2, if both of the functions \f/ f and \p" satisfy 
the boundary condition of (1.4) on these surfaces. 

We here introduce the solution \p m (x, x/ ) of (1.3) for the smooth and homogeneous earth 
having the propagation constant k m . By the use of (1.5) and also the boundary condition, 
the symmetrical relation of yp, t ,(x, x f ) is readily derived; putting \p* (x)=\(/ m (x ) x x ) and \f/"(x) 
= \fs m {x,x 2 ) in (1.5), the points x { and x 2 being arbitrary points in space, and taking the whole 
space for 2, we have 

^m(«l, X 2 )—yp m (x 2y X 1 ) = — [\fs m (s, Xi), lpm(s, X 2 )] Sm =0, 

or 

yp m {xu x 2 )=\p m (x 2) ./;,). (1.7) 

It is appropriate to begin with the simple terrain as in figure 2. Here, the propagation 
constants k 2 and & 3 of the ground are assumed to be different. 

Now in the preceding Green theorem (1.5), let \p / (x) = \f/ 3 2(x^ x) and ^ ff (x)=\l/ 2 (x ) x } ). Here 
yf/ Z2 is the Green function to be obtained, and \f/ 2 is the Green function for the homogeneous ground 
of the elevation z=a 2 and of the propagation constant k 2 . Taking the space 2 as the whole 
space except the space enclosed by the surface s 3 defined in figure 2, and the point x x being in 2, 
we have (a 2 >a 3 ) 



or 



fefe, x { ) — ^Ax u ./J + ifefe, s z ), \l/ 2 (s 3) Xi)] t z±>a 2 

^82 (»4, a5l) = [^32(^4, S3), lfe(*8, &l)]> Zl<a-2, 



(1.8) 
(1.9) 



when the point # 4 is below the surface s 3 . Here, it is noted that the integrand of the surface 
integral is continuous at the boundary of the medium k 2 and the atmosphere because both 
1^32 and \p 2 should satisfy the same boundary condition, and therefore there is no contribution to 
the surface integral on the right side of eqs (1.8) and (1.9). 

Just in the same way, by exchanging the roles of the Green functions \[/ 2 and ^ 3 , we have 

(a 2 >a 3 ) 

fe(x 4 , xO=\p s (x h a?i) + [^82(a?4, 82), h(s 2j Xi)]. (1.10) 

These two equations are complementary to each other, and, by the method of successive 
substitution, we have the required solution in series 



fe0&4, x l )=\l/ 2 (x 4 , ajO + hfofej s z ), t// 2 (s 3j Xi)] 

+ \MX4, 82), fofe, *s), fa(s Z , Xi)] + 



4 >a 2 . (1.11) 



Figure 2. The form of terrain and the integration 
domain for eqs (1.8) and (1.9). 




41 



I I(M<\ 



[+, v, *"]=[[*, n +"]=[+, w,+"]], 



and, on the vertical boundary surface, the boundary condition (1.4) is taken into account. 
On the other hand, when the points x { and x± are located over the different earth media, 
say on the sides of the media k 2 and k 3 , respectively, the term fefe, S3) in the bracketed terms of 
(1.8) and (1.9) may be given by 

fefe, s 3 ) ^ 3 to, 5 3), 

except in the vertical domain of s 3 which is not important for the integral, because \// 2 (ss, £1) 
rapidly tends to vanish with the distance a 2 — z by the term exp \—i\k\—k\{a 2 —z)\ in the 
ground. 
Thus 

fete, ffi)~ife(*4, frO + W'sfe, 5 3 ), fa(ss, xi)], Z4>a 2 

^#3(^4, s 3 ), hiss, &i)L Zi<a 2 . 

The above result is just the equation obtained from (1.11) by the omission of terms of 
higher order than the third in the series. The higher order terms become important only when 
the point x A or x x is located in the immediate vicinity of the vertical boundar3 T surface. Also 
the surface integral over the vertical boundary surface can be neglected because of the large 
attenuation of ^ 2 fe, #1) in the ground, and so far the presupposition that the horizontally 
polarized wave will not be induced is self-consistent. 

There may be another problem of singularity of solution at the diffracting edge or unique- 
ness of solution for the form of boundary assumed. But it is known [Born and Wolf, 1959] 
that, insofar as the field components can be expressed as convergent Fourier integrals, which 
is the situation in our case, they are free from having singularities of too high an order, and their 
uniqueness is assured. 

The above result can easily be extended to the general case where there are several bound- 
aries of the earth medium discontinuities, as in figure 3. As is illustrated there, we successively 
assume, along the wave path from the point x x to x n+1 , the values k 2 , k z , k±, . . ., k n for the 
propagation constants of the different sections of inhomogeneous earth, and the values a 2 , 
a 3 , a 4 , . . ., a n for the radii of the surfaces of the respective sections, and r 2 , r 3 , r 4 , . . ., r n for 
the propagation distances, respectively. Also, the solution in this case will be denoted by 

\p n 2 (x n+ i, x x ) and thus the solution \f/ n _i 2(^+1, #1) will be the one to be obtained from 

fn 2(^+1, Xi) by setting k n ->k n ^ and a n -^a n . x . 



On referring to figure 3, we now set \f/ f (x)=\f/ n) 



2\X n 



+"(x)=+n- 



2 (a?, xO in 



eq (1.5), and take, for the space 2, the whole space excluding the space of the medium k n 
and the atmosphere above it in the range z<jL n -u the point x x being assumed to be within 2. 
Then, as in eq (1.8), we have 

$*..... 2(Xn+l,%l)=tn-l,..., 2^71+1 , «l) + [fn 2(^+1, «n), tn-1,..., 2(«», &l)], 3»+l>«n-l 



or 



(1.12) 

tn,...,2(%n+i, #i) = ty'*,...,2(a;»+i, s n ) , ^„_ lf . . . t 2 (s n ,x 1 ) ], z n+1 <a n - U (1.13) 

depending on whether the point x n+i is within or without the space 2, respectively. 




Figure 3. The form of terrain and the domains of 
integrations for eqs {1.12), (1.13), and (1.14)» 



42 



Id the same way, by letting ^ r {x)=yf/ n 2(#»+i, x) and ^"(#)=^»(a;, Xi) in eq (1.5), we 

have 

^n,..„ sC^iM-li »l)=^»(»iH-li ^l) + [^n 2(^+1, *'), **(*', »l)l- (1-14) 

Here §' is the whole upper surface of the media k rn for which a m >a n , m = 2, 3, . . ., n— 1, 
plus the surface of z=a n above the media for which a m <^d n , plus the vertical boundary surface 
of the medium k n , as is illustrated in figure 3, and both the points X\ and x n+l are to be outside 
the surface s'. 

When the solution ^ n _i 2(x n+1 , #i) is known, we can, by the successive substitution 

of eqs (1.12) and (1.14), express the solution f„ 2(^+1, %i) in a form of series similar to 

eq(l.ll): 

^n 2(^*1+1, »l)=^n-l f ... i 2(^»+l, Xl) + [^n(Xn+l, S n ), ty n -\ 2(*»> ^l)] 

+ [**-! 2^+1, O, *»(«', Sn), *-l 2(*.,*l)]+ • • ., (1.15) 

or, in just the same way, 

^,...,2(^+1, 3i)=[^»(z»+i,«»), l^ii-l 2(s n ,x 1 )] 

+ [A.(*»+l,«.),*.-i.....s(^^ . . -, (1.16) 

depending on whether the point a? n+] is within or without the space 2, respectively. Thus, 

eq (1.15) gives a recurrence formula for \p n 2 . 

As in cq (1.11), we can neglect the terms of higher order than the third or second in the 
series above except for tin 4 range immediately near the boundary. 

2. Formula of Field Strength in the General Case 
The attenuation coefficient A is here defined by 

$n 2(^+1, x 1 ) = 2A(z n+l \r /n r„_ l r 2 \z 1 )^ (x N + u x x ) 7 

Mxn+u ^=11} e ~ ikir > r=r n +r n ^+ . . . +r 2 , (2.1) 

where \l/ (x n+u x } ) may he regarded as the solution in free space. Then, on the assumption 
that the boundary surfaces between the different sections of the in homogeneous earth are 
all vertical and parallel with each other, the result of evaluation according to eq (1.15) is as 
follows (figure 3): 

A(z n+1 \r H ,r n - ltmmmt r 2 |zi)= S ^lr n ) h A{z n+l \r n ) t T{r n . l ) tnttn _ l 

XT(r fl _^ tti _ vtn _ 2 . . . nr z ) t ^ z T{r 2 ) H , t J t2 (z u a 2 ). (2.2) 

Here, t m stands for the set of roots of the equation 

W'(t)-q m W(t)=0, (2.3) 

where 

PF(-*)=exp (-i2*/3)(irt/S)»ff#(§ t 3 '^ 



fk^kl-mi, Vert. Pol. 

*ft»=(M2)»X-j 

LV*t— *S/*i, Horiz. Pol. (2.4) 

These notations are similar to tliose of Fork [1946]. 

The factor /, 2 (zi, Ch) is the ordinarily defined height-gain function for the end point x t : 

f t2 (zu<h)=W{h- yi )IW{t 2 ), yi=(2/k 1 a)ik 1 (z 1 -a 2 ), (2.5) 
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and the other two kinds of terms that appear in the preceding formula for attenuation co- 
efficient A are as follows: 



4(2n+ik„)j n =V(ir/2)fc^^^ 

(2.6) 

f 2m/ J m («« , dm) — <ln ft m (d n , O , ^n > d m 

Xexp [-i(rja) {h(d n -d) + (k 1 a/2)H n }]X < (2.7) 

t 2m/*,, (&m , dn) ~ gn/ f B (^w > <0 > ^m > <&* . 

The term yl(,3 w+1 |r„)^ depends only on the electrical properties of the nth section through 
q n and t n and on the propagation distance r n and the elevation a n of the nth section, but does 
not depend on those quantities of other sections. 

On the other hand, the term T(r n ) t >t depends on both the mth and nth sections through 
t m and t n and also through q m and q n , and therefore it serves as the coupling term between the 
two sections. Besides, it depends of course on the height difference between the two sections 
through the ordinary height-gain function f t (a m ,d n ) and also through another height-gain 

function f t (a m , d n ) defined bv 

j i n \ in) n j j 

Uftjdn, d m ) = — (hd/2)* „ ftjfln, d m ), 



dfcaj- 



(2.5 



which is proportional to the first-order derivative of the ordinary height-gain function. 

There is a clear one-to-one correspondence between the terms of the formula (2.2) and 
the respective sections of the inhomogeneous earth, i.e., the height-gain function ft 2 (zi,d 2 ) 
expresses the effect of the elevation of the point x Y from the ground, of course, and the term 
T{v2)hm expresses the effect of propagation along the surface of the section of No. 2. In 
the same way, the term T(r z )t it t 3 corresponds to the propagation along the section of No. 3, 
the term T(r^) h ,u corresponds to the propagation along the section of No. 4, and so on, and 
finally, the term A(z n+l \r n ) tn corresponds to the propagation along the nth section, over 
which the point x n+1 exists. 

The convergence of the multiple series of the attenuation coefficient A is ver}^ good 
when the propagation distance of every section of the inhomogeneous ground is sufficiently 
long, and in this case, the first term of the residue series is a sufficiently good approximation. 
This situation is just the same as in the ordinary Van der Pol and Bremmer formula for a 
homogeneous spherical earth. 

On the other hand, the convergence of the series becomes poor when the propagation 
distance over one of the sections is very small. In this case, the flat-earth approximation or 
other proper approximation can be suitably used. 

In the extreme case where the width of the one of the sections tends to zero, it will represent 
a ridge on the ground, as in figure 5. In this case, the responsible series in the formula takes 
the asymptotic form as r m tends to zero (fig. 4), 

lim ^t m T{r m ) Ht t m T(r n ) tmt tfl = T (m) (r n ) H , tn = {qf tl (d m , a)j tn (a mj a n ) - q n j h {d m} a t )f tn (a m9 a n ) } 

r m —>0 

X{k l (d l —a n )(2/k 1 a)^+t l —t n }~ 1 (t n —ql)- 1 exp [—i(r n /a){k 1 (a n —a) + (k 1 a/2)H n }], (a m >a nj a). 

(2.9) 




Figure 4. The form of terrain and the notations for 
eq {2.9). 
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The above result is applicable to the case of diffraction by ridges on spherical earth, as 
it will be seen in the example of the following section. 

It may be remarked that, in the special case of a n =a n ^ l = . . . =a 2 , the formula (2.2) 
reduces to that previously obtained for mixed paths on a smooth earth [Furutsu, 1955]. 

3. A Few Examples of Applications 

The simplest application of the foregoing formula (2.2) will be to the case of one ridge on 
a spherical ground, where the electrical properties on each side of the ridge could be different, 
as illustrated in figure 5. By the use of the result (2.9), the attenuation coefficient A is then 
given by 

A=^ UJ ^lr^A{z h \r^ u T^{r 2 ) Kh j t2 {z u a 2 ). (3.1) 

One of the immediate applications of this formula will be the estimation of the effect 
of a ridge on a spherical smooth earth (fig. 6), and another will be the similar estimation of the 
effect of a cliff on the ground radio wave (fig. 7). In the case in which the propagation distances 
/' 4 and r 2 on both sides of the ridge (or cliff) are long enough and the height h of the ridge (or 
cliff) is low enough, i.e., 

&i(r 4 — \ 2ah) {k x a)-^ 1 , /•,(/•,- - N 2ah) (£,a)^» 1 , 

so that simply the first term of the series would give a good approximation, then the effect of 
the ridge (or cliff) can be expressed by multiplying the field strength in the case of smooth 
earth by the factor Fif. Here, using the notation t" and t" for the first values of the set of 
values t A and t 2 , 



n 



<h /'/<.(« 3 , «.4)./;.>(«3, a 2 ) -fc/ ( o(a 3 , a 4 )/;»(a 3 , a 2 ) } X {^(a 4 -a 2 ) (2/*,a)*+tf-/2} -'(tt-gf)- 1 (3.2) 



or, when a 4 =a 2 and (ji = q-2, 



(3.3) 



Tims, depending' on whether the boundary is a ridge or a cliff, the field strength is given 
in the form (figs. 6 and 7) 



K= {E) h , F KK , (p) , E = (E) h=0 F$ K , ( P ) , 



(3.4) 



respectively. Here, K and K' correspond to k and ¥ ', respectively, and according to K. A. 
Norton, K and b are defined by 



and 



K= \q\~ l 2- ll \ -b = r/2-\-2 arg (< £ ), 
p = k l h(k 1 a)~ 1/3 , ki =ca/c. 



(3.5) 
(3.6) 



Hence p is proportional to the height of the ridge or cliff. 

The set of graphs in figure 8 shows the numerical values of the factor F KK , giving the effect 
of a ridge on the ground waves, and it is displayed as a function of p, which is proportional to 
the height of the ridge, for a useful range of values of K and b. 



Figure 5. The form of terrain and the notations for 
eq (3.1). 
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Figure 6. The form of terrain for eq (3.3) and 
figure 8. 



Figure 7. The form of terrain for figure 9. 
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Figure 8. Ridge effect on ground radio wave as a function of p and b in the case of K'=K. 



In the same way, the set of graphs in figure 9 shows the numerical values of the factor 
Fkk giving the effect of a cliff on the ground wave, and it is also displayed here as a function 
of p for the same range of K and b. 
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Figure 9. Effect of cliff on ground radio wave as a 
function of p and b in the case of K'=K. 
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4. A Few Examples in Flat-Earth Approximation 

The results of eq (8.4) are not available in the case where the propagation distance on 
one or both sides of the ridge is very short. Then the effect of the ridge depends on the propa- 
gation distances from the ridge, so the result is not so simple as (he preceding cases. 

In this case, the effect of the earth curvature may be neglected and the flat earth ap- 
proximation is appropriate; the convergence of series of the formula (3.1) becomes very poor 
in this case. However, there is a general rule to overcome this difficulty which will be treated 
in part II. The results in the flat-earth approximation are then obtained as the asymptotic 
forms for a-^co ? i.e., the infinite radius of earth. The formula (3.1) in the case of a flat earth 
becomes as follows: 

Referring to figure 10, it is convenient first to introduce the notations that will be used 



later: 






/4=(a 3 -a 4 )v / *i/2/'4^ ,r/4 , 






a 2 ) V^i/2r 2 e /7r/ S 
(a z -a 2 )-y/kJ2^ 2 e iv/ \ 



h 



oo 7 , I kj-yjkn k 1 

— —; k n = < 

c l k\l^kl-k\ 



klNk\-k\, Vert, Pol. 



Horiz. Pol. 



(»=2,4). 



(4.1) 



Figure 10. The form of terrain and the notations for 
eqs (4.2) and (44)- 
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The attenuation coefficient in this case will conveniently be expressed by A (3) (z 5 \r h r 2 \zi), 
which means the attenuation for the wave propagation from the point x x at the height z x to the 
point x 5 at the height z 5 across the ridge of the elevation a z which exists at the distance r 2 from 
the point X\ and the distance r 4 from the point x 5 . 

In the special case where both the points x x and x- are on the ground, i.e., when s 5 = a 4 , 
z x =a 2 orj 5 =f 1 =0 } 

X ydid 2 /n 4 n 2 £ } (f4+ i VdJ S l (f 2 +i V^ 2 ) + (djn^{& (^4X^2/4+ V^j^M) 
~<f (p 4 , (/ 4 +iV^)/p 4 )}+ feM 2 ){#((r 2 , Un 2 f,+^ A ;f 2 )l<j 2 )-S( P2 , (J2+1JQIP2)}]]. (4.2) 
Here 

rf (z,n/z)=(^)\ z2+n2 ( m dxe~ xi P <T^, 

W^V Jz J (n/z)x 



\ 



W/4— y/fyfa+i Tjdjrii, p 4 =i ^'d i n 2 /n 4 —f 2 , 



a 2 = \n 2 f 2 — in^f i -+i^d 2 /n 2 , p 2 = i\d 2 njn 2 —f 4 . (4.3) 

In the general case of z$^a 4 and z x ^a 2 , 

A^\z,\n,r 2 \z 1 )=F(d i J i +f 5 \d 2 J 2 +J0 
-^(<k/d-/ 5 |0J a 4^ 

(4.4) 



Immediate applications of the formula (4.2) are possible for the effect of a ridge on the 
ground wave, and the effect of a cliff and the effect of a bluff along a coast line, as illustrated 
in figures 11 to 13, some numerical results of which will be displayed in the following. The 
notations that will be used are: 

h : Height of ridge, bluff, or cliff, in meters S 

/ : Frequency in megacycles per second 

a : Conductivity of land in millimho per meter 

£}=!?x{:;' /<i*/.. (4.5) 



Case 1. One Ridge on Homogeneous Ground 



The ridge is assumed to have the height h from the flat ground and the transmitting and 
receiving points to be at the distance r 2 and r x from the ridge on opposite sides, as illustrated in 
figure 11. 

The two curves in figures 14 (a), (b) represent respectively the attenuation coefficient 
and the phase delay in the case of a homogeneous earth. Here, the abscissa is the Sommerfeld 
numerical distance for the whole propagation distance. On the other hand, the curves in 



Figure 11. The form of ridge and the notations for 
figures 15 (a) and (b). 
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figures 15 (a), (b) represent the same attenuation coefficient and the phase delay in the case 
where one ridge is present between the transmitter and receiver. Here, the parameter d is the 
numerical distance on one side of the ridge, and ) r is the numerical height introduced in eq 
(4.5) which is proportional to the height of the ridge. 

Case 2. A Bluff at Coast Line 

The bluff is assumed to have the height h along a coast line, and the transmitter and re- 
ceiver to be at the distance r\ on the sea side and r 2 on the land side from the coast line, as 
illustrated in figure 12. The sea is assumed to be a perfectly conducting plane, while the land 
has the finite conductivity a. 

The curves in figures 16 (a), (b) represent the values of the attenuation coefficient and the 
phase delay in the case where there is no bluff at all at a coast line. The abscissa S is the Som- 
m erf eld numerical distance for the propagation distance on the sea side measured in land con- 
ductivity, and the parameter L is the numerical distance on the land side. 

On the other hand, figures 17 (a), (b) show the values in the case where there is a bluff 
at the coast line. Here Y is again the numerical height proportional to the height of the bluff. 

Case 3. A Cliff 

The cliff is assumed to have the height h on the flat ground, and the transmitter and re- 
ceiver to be at the distance r x on the lower side and r 2 on the higher side, as illustrated in 
figure 13. 

The curves in figure 18 display the corresponding values in this case. 

In case 1, ridge diffraction, the magnitude A of the attenuation coefficient decreases with 
the height h of the ridge in most cases, but for the large transmission distances of D^>1 it 
ceases to decrease at some height and then tends to gradually increase witli the height. This 
fact may be interpreted by noticing that the diffraction loss by the ridge is rather smaller than 
the transmission loss along the dissipative ground at large transmission distances; a ridge could 
give an obstacle-gain even on a flat earth. 

In case 2, where the radio waves propagate 1 across a coast line having a bluff, the rate of 
change of the relative phase with the distance r { (which is proportional to the numerical distance 
Shi figure 17) from the coast line becomes larger as the height of the bluff increases and is some- 
times much more than that without the bluff. 

Besides the subjects just mentioned, it is also possible to evaluate the diffracted wave 
by a ridge of finite thickness, as illustrated in figure 19. 

In summary of part I, it may be remarked that the general formula (2.2) for the attenu- 
ation coefficient, given in the residue series, takes the unified form independently of the number 
of sections of the inhomogeneous earth, but the corresponding formula in the flat earth approxi- 
mation is probably more difficult to derive and would take a more complicated form. Equation 
(4.2) or (4.4) is an example of the latter for two sections. Practically, the case of three sections 
will be the limitation in which the flat earth approximation is possible, unless some special 
assumptions about the propagation distances are made. Also, the effect of reflecting waves 
from the boundaries of discontinuities on the transmitted waves are completely neglected 
here. But they are generally believed to be very small in most cases. In fact, the result 
has been proven to be exact at least for the trapped wave (or surface wave) mode [Furutsu, 
1959b]. 
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Figure 12. The form of coast line and the notations Figure 13. The form of cliff and the notations for 

for figures 17 (a) and (b). figure 18. 
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Figure 14. The attenuation coefficient (a) awrf the 
phase delay (b) in the case of a homogeneous earth. 
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Figure 15. An example of the attenuation coefficient (a) and the phase delay (b) in the case of one ridge present 

Part II. Multiple Diffraction of Electromagnetic Wave by Surfaces Having 

Finite Radii of Curvature 

1. Statement of Problem and the Result in the Simplest Case of One Diffracting Surface 

Many authors have discussed the diffraction of electromagnetic waves by a single object 
such as sphere, cylinder, paraboloid, ellipsoid, plane, etc. However, there are many difficulties 
when trying to calculate the diffraction loss by actual mountains or hills, because they usually 
can not be seen as a single smooth diffracting obstacle, but more or less have irregularities. 
Strictly speaking, the actual diffraction of radio waves by mountains is a problem in multiple 
diffraction. But it is, of course, impossible to derive the formula of diffraction loss for such 
complicated obstacles, and different models have usually been adopted for the diffracting 
mountains according to the different range of frequency. However, if we could assume that 
every diffracting mountain has a smooth surface whose radius of curvature is sufficiently 
large as compared with the wavelength, the problem can be treated on fairly general conditions, 
and we can get the result in a unified form, independent of the number of diffracting mountains. 
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Figure 16. An examvle of the attenuation coefficient (a) and the phase delay (b) in the case of no bluff at coast line. 
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Figure 17. An example of the corresponding attenuation coefficient (a) and the phase delay (b) in the case of a 

bluff present. 



In the simplest case of one diffracting surface, the result can be deduced from the ordinary 
Van der Pol and Bremmer formula for diffraction by a large spherical surface. Referring to 
figure 20, the attenuation coefficient A, which is the ratio of the field strength to that in free 
space, is given by 



A=^kX2/k 1 ay^(t s -( 1 2 )- 1 ex V [-i{(r/a)^ (1.1) 
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Figure 18. An examvle of the attenuation coefficient 
(a) awrf tAe p/mse deZa?/ (b) w */ie case 0/ one cZz^ 
present. 
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Figdbe 19. The form of ridqe having finite thickness. 



Figure 20. The spherical diffracting surface and the 
notations for eq (1.1). 




Here (,'s are the roots of 

W'(t)-qW(t)=0, W(-t) = ^e- t2r/ Hm% ) (it*), 

and 

f.(y)=W(.t,-y)IW(U), y i ={2lk l a)^k l h u %<a) 

(k^'kf^Fjkl Vert. Pol. 
Pol. 



\ k^'kt-mi Vert. 
''2=(W2)«X , . 

i yki—Ki/ki, Honz. 



For 2/^1, the height gain function /,.(?/) takes the asymptotic form 
f s (y)^y-y* exp [-iiiy^-tsy^+iy-ynl+irm/Wits), 



(1.2) 



(1.3) 



(1.4) 



where the argument of the exponential function is expanded with respect to t s and the terms 
involving powers higher than 2 are neglected. Thus, when the heights of either or both the 
transmitter and receiver take sufficiently large values as y u y£$>l } the attenuation coefficient 
A of (1.1) takes the following forms: 
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Case 1: ?/i»l,y.><l 
/l~exp [- i (ha/S) (n/ay\2^ A V/n 2 (*. - <f ) ~ 'W ' (f .) " ' X exp [-i { 0(^/2) w*. +^/2 } ]/, (j/ 2 ) , 

d=(r—r 1 )/a, 7\=^2ahi ) r 2 =^2ah 2 . (1.5) 



Case 2: yi»l, // 2 >1 



—[-¥{©•+©"}] V^p* 



(kia/2) 1 '* (t t -f)-W(fi 9 )-* 



Xexp[-t{0(£ 1 a/2) 1 / 3 * s +(2^^ 

e=(r—r 1 —r 2 )/a. (1.6) 

In case 2, the series is convergent for the whole range of 0. On the other hand, in case 1, 
the series diverges for negative values of 0, and hence it must be analytically continued, as will 
be shown later. Thus the analytically continued function A of d is regular on the whole range 
of 0, including negative values. However, in case 2, if the square term of t] were neglected in 
the exponential function, the 4 function A of would have a pole at 6 = 0, even though it is 
analytically continued. Thus the square term of f; cannot be neglected in case 2. 

The convergence of series of these formulas is good for large diffraction angles. But, 
otherwise, it becomes poor. However, there is some general rule to overcome this poor con- 
vergence, as will be treated in the following section. 

Since the above results have been developed according to the original Van der Pol and 
Bremmer formula given in residue series, they could he valid only when both the transmitter 
and receiver are not far from the surface, and it is not immediately clear whether they arc 
also valid even when calculating the field far from the surface. 

But it can be proven [Furutsu, L956; Wait and (onda, 1950] that the results (1.5) and 
(1.6) are correct even when the transmitter and receiver are at great distances from the dif- 
fraction surface as compared with the radius of curvature, if we reinterpret i'i(^^2ahi) and 
t 2 {t £ -yJ2ah 2 ) as the lengths of the parts of wave path from the transmitter and receiver to the 
first contacting points on the diffracting surface, respectively, as illustrated in figure 20. The 
important part of the mountain surface which decisively contributes to the diffracting waves 
is the very small part of the surface in the vicinity of wave path, and the other part is not 
important, provided that the radius of curvature of surface is sufficiently large as compared 
with the wavelength. 

On the other hand, in this small part of the diffracting surface, the surface 4 could be ex- 
pressed as a surface of second degree having some finite radius of curvature such as a sphere, 
cylinder, paraboloid, etc., and therefore it follows that we could solve the wave equation 
exactly in the range of the important part of the diffracting surface. In fact, the Green func- 
tion usually happens to be the same in this small range, independently of the kinds of surfaces 
adopted, provided that the radius of curvature of the diffracting surfaces is defined along the 
wave path. 

2. Methods for the Case of Poor Convergence 

For kinds of series such as (1.5) and (1.6), there are some general rules to overcome the 
difficulty when the convergence of series is poor. Taking into account 

we have for the arbitrary series 2 s a(£ s ) 
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where (7 S is the infinitesimal contour integration path around the sth pole. 

Hence, if the integrand thus formulated does not have any pole besides the poles t s 's, as 
in case 1 of eq (1.5), the sum of the contour paths is equivalent to the contour path C around 
the set of poles, as illustrated in figure 21. Further, if the integrand tends to zero sufficiently 
rapidly at infinity, we can deform it to the path Ci~\-C 2 in figure 21; these paths are proven 
to be the best paths for numerical integration in the meaning that the integrand decreases 
most rapidly on these paths. Furthermore, in the case of (1.5) the integral converges for 
negative values of 0, even when the original series diverged. Thus 

i:f = f=f ■ (2.4) 

s JC S JC JCi+02 

On the other hand, however, if the integrand has extra poles, such as those of the func- 
tion W(t), besides the necessary poles of t s , the foregoing method cannot be used, and another 
integrand must be sought. 

From the Wronskian identity, we see that 



]qv(t s )-v'(Q}W(t s )-- 



=1, v{-t)=\ Jj {e-'WHtft Q t^+e^H^mh (f * 3/2 )}« (2.5) 



Hence, in principle, we could multiply the integrand by this function to any power we wished. 
For example, by multiplying it once, eq (2.2) is replaced by 

? a W =S ?X, ( w'U-mti ) {t -^(t)am. (2.6) 

Thus, even if the old integrand had the undesirable extra poles of the function W(t), 
the new integrand would not have them. Thus we have the contour integration path C and, 
further, the path d + C^ if the integrand decreases sufficiently rapidly at infinity. This 
situation actually occurs in the series of (1.6). 

Here the question may occur whether it is possible to multiply the integrand by the 
square of the function (2.5) and to deform the integration path to the path Ci + d. The 
answer is no, because, though the integrand has no extra pole besides t s , it diverges at infinity 
on the way deforming the path from the C to the path Ci~\-C 2 . Hence, generally, there exists 
only one integrand for which the integration path Ci-\-C 2 is available. 

Usually the Kirchhoff approximation terms appear as the leading terms of these integrals, 
and are obtained from the asymptotic forms of the integrands for large magnitudes of t; for 
instance, 

-1 f v'(t)-qv(t) \ _i ( 

2wi \W'(t)-qW(t) J 4tt ^' n 

on the path d or any path of (0, oog -i 0) in the range x>|8>7r/3, and it tends to zero as |£[— >°? 
on the path (\. 



Figure 21. The infinitesimal contour path C s and 
the integration vaths Ci and C2 for eqs {2.2) and 

CM). 
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Using these asymptotic forms in the integrand of A for case 2, it becomes 



yl^exp [ ] 



Jr-J- (ha/2)^ -L P '\lt exp [-i{^(A: 1 a/2) 1 ^ + (2i 1 )- 1 (l/r 1 +l/r 2 )(A: 1 a/2)2/3^} ] 
=exp[] i=^ 2 r°° , *" |2 <fc, f=^V^ir 2 /2(r 1 +r 2 ), (2.8) 



with 



exp []=exp [-i(k 1 a/3){(r 1 /ay+(r i /ay}]. 

The result of exact evaluation takes the form 



A=exp [] \^e H ^ )2 f e- t2 dt—nG{i)\> 



$=6(hal2yi\ ^(ha/2)"* V(2/*i)(l/r 1 +l/r 2 ). (2.9) 

Here, the numerical values of the function 6?(£) have been calculated by Logan [1959] 
for a surface of perfect conductor, and by Wait and Conda [1959] for a wide range of surface 
impedance. Tt may be remarked that the corresponding numerical integration method for 
the case of poor convergence has been used more previously by Fock [1946] and Rice [1954]. 

3. General Formula for Diffraction by Two or More Surfaces 

As illustrated in figures 22 and 2:>, we now consider the general case where the waves 
propagate over several mountains with the subscripts 2, 3, ...,// along the wave path 
from the point x { to x ti+1 , and the radii of curvature (along the wave path) a 2 , a 3 , . . ., 
a n , and the propagation constants k 2 , k z , . . ., k n , respectively. The diffraction angles 
of the respective mountains will be denoted by 6 2 , 6 :h . . ., 0,„ and thus d m =a m $ m (m=2j 
3, . . ., n) will be the distance of that part of the wave path contacting the ?nth mountain. 

Now, on the assumption that the whole wave path lies in a plane profile, the result of 
evaluation gives the following expressions for the attenuation coefficient A which are just 
the generalizations of (1.5) and (1.6). Using the notation t m for the set of roots of eq (1.2) 
in the case of q=q m , they are [Furutsu, 1956]: 

Case 1: Wr m ,, ± i) 3/4 (iiO" 1/4 «l (2<m<n) 



A — i Vn+l, n~\d n -\- 



+r Z2 +d 2 +r 21 )/k? l r n +i tn r n>n -i . . . r Z2 r 2] 



X 2—i *■ Vn+l,n)o, *--*(£») t n 

t n h,h 



m,) t J(r, 2 ),'m 2 ) t T(r 2l ) t , 2 , {) . (3.1) 



Case 2: (a 2 /r 21 ) 3/ \k ] a 2 )- 1/4 > 1, (a 2 /r 32 ) 3/4 (it 1 a 2 )- 1/4 <l, 
(amAVm±i) 3/4 (#iO~ 1/4 <l (3<m<ri) 

A={(r n + hn +d n + . . . +r Z2 +d 2 ) /kr 2 r n + lin r n , n -i • • • r AZ r z2 } m 

X S T(r n+ltn ) ,t n TUn)t n . . . T(r^ t4 , H T(^ tz T(r n ) H , t$ T(b A ) h . (3.2) 

tn • • ., H, tl 





Figure 22. The form of diffracting surfaces and the Figure 23. The form of terrain and the notations for 

notations for eq (3.1). eq (3.2). 

55 



1/2 



Case 3: (<W*i) , ' 4 (*i<h)- i ' 1 >1 j (a„/r K+1 , B ) 3 ' 4 (^ B )-"*>l 
(«„/r ra , m±1 ) 3/4 (^a m )- 1/4 «:l (3<m <»— 1) 

o4={(<4+r B , n _ 1 + . . . +rf 3 +r 3 2+rf2)/^r 3 ? , „, n _ir„_j,„_ 2 . . . r i3 r 32 } 

X 2 T(l».3 B+ i)« B T(r B ,„_ 1 ) (B , % _ 1 T(?„_ 1 ) (B _ 1 . . . T(^ H T(r^ H , h T(b.edtr (3.3) 

Here, r m>n is the distance between the mth and ?ith mountains and 

km — Vm \k\dml2) , 

T(r) tm , ln =exp [-i(2k l r)- 1 {(k 1 aj2) l »t m -(!c 1 aj2) 1 >H n }*\, 

TU, z) tn =2Mt m -<?J- 1 W(t m )-y t Jz, a,,)e-^'^n>. (3.4) 

In case 1, both the transmitter and receiver are sufficiently apart from the diffrac ting- 
mountains (fig. 22). The formula (3.1) in this case consists of two kinds of terms, one of 
which, T(l- m ) t , depends only on the diffraction angle, the surface impedance and the radius 
of curvature of the mth mountain, but does not depend on those quantities of other diffracting 
mountains. 

On the other hand, the other kind of term T(r mttl ) t ,t n depends on the radii of curvature 
of the mth and nth mountains and their electrical properties through t m and /„ and also the 
propagation distance between them, but does not depend on other quantities. Therefore, it 
serves as the coupling term between the mth and nth mountains. There is no other kind 
of term. 

This fact facilitates the understanding of this formula considerably; there is a clear one- 
to-one correspondence between the terms in the formula and the respective parts of the wave 
path: for instance, T(r 2 i)t 2 ,o corresponds to the wave path from the point Xi to the second 
diffracting surface, the term T(£ 2 )t 2 corresponds to the wave path along the same diffracting 
surface, and T{r 32 )t s ,t 2 corresponds to the propagation in free space from the second to the 
third diffracting surface, etc. 

As already stated, this formula is essentially the generalization of the formula (1.6) for 
one diffracting surface which was derived from the original Van der Pol and Bremmer formula 
on the restrictive condition that both the transmitter and the receiver are not so far from the 
diffracting surface as compared with the radius of curvature. But here it is derived from this 
general formula as the special case of one diffracting surface, on the more general condition. 

In case 2, one of either transmitter or receiver, say the point x u is on or near the diffracting 
surface (fig. 23). The only difference from that of case 1 is in the last two terms; they were 
replaced here by the term T(£ 2 , Zi)t 2 > and the others are just the same. Here, ft m (z, dm) is the 
ordinarily defined height-gain function. Again, the special case of this formula agrees with (1.5) . 

In case 3, both the transmitter and receiver are on or near the respective diffrac ting- 
surfaces. The change to be made in the preceding formulas is so evident that it might not 
be necessary to mention it here explicitly. 

Sometimes it may be necessary to take into account the contributions of the waves which 
are reflected from the surfaces between the diffracting mountains. Mathematically speaking, 
these reflecting points on the surfaces are just the phase stationary points of the integrand, 
and we could get the result from these formulas by multiplying the reflection coefficients of 
the surfaces and also adjusting the diffraction angles so that they correspond to the reflected 
wave paths. 

4. Summation of the Series 

Finally, there remains the problem of convergence of the formulas (3.1) to (3.3). There 
are two kinds of series in the above formulas; one is of the form 

s^ii.j^+^i/yij'Ovm-i)/,^^!, (4.i) 
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and its convergence is poor when £,„ is near or smaller than 1. Now the treatment for the case 
of poor convergence is exactly the same as in section 2. Using the same method, we have 
the result 

eq (4.1) = 2 r? (r w+ i. m)t m+1 ,oT(r mim . 1 ) ,t m _ l Mha m /2)^F^ m -^ rj m , q m ). (4.2) 



Here, on referring to eq (2.9), 



Vm= 0W2) 1/3 V(2/&i) (l/r M+1 , w +l/r m , w _0, 

^(fca^'Mft^^ (4.3) 

and % m depends on the preceding and succeeding diffracting surfaces through t m ±u and thus 
it serves as the only coupling variable between the diffracting surfaces. 

When the diffraction angle of the wth surface is small and/or the radius of curvature is 
sufficiently small, the second correction term becomes very small in most cases and can be 
neglected as compared with the first leading term; in the case of£ TO , rj m <^] the series (4.1) takes 
the simple form 

S^ 1 ^ ' • m)t m . .oT^.Oo, ^. 1 W({ w -{:)/iI w «'' /4 )Vfc/(l/r» + U+l/^»-l), fm<l , 

£(z) = ^.e' 2 ( m e- t *dt. (4.3a) 

Another special case is the case in which the diffraction angles of the m±lth mountains 
are sufficiently large, as f m ±i/i7 m ±i^> 1 . It follows then that £,Jr} m <^\ and thus £'„ in (4.2) 
can be neglected: 

]E=r(r w+li J lm+li or(r^^ (4.3b) 

Similarly, on the condition of f w /^ w »l and £jrii> — l(Z^m), 
^(Jfcia^) 1 ^(^-^,i ?lllf2 J^V5(Jfc 1 aJ2) 1/8 \- I -e-^-GiU^MiUarn), W>Kl (4.3c) 

L 2 \ *£« J 

Here, the function M(£ m , a m ) is introduced for later convenience. 
The other kind of series is 

S7'()- ;!2 ), s ., 2 Tfe, *i) ,2^r(r 32 ) h , oS^fe-f,, gt)h, (4.4) 

which always occurs with the height-gain function of receiver or transmitter. Using the same 
method as in section 2, it takes the form, on the condition that £i/r}i^>l(lj*2) and &/%>— 1, 

Sr(f2-«,2,) h ^sf(&;{i+i(«/*i)2i}, (2/fca 2 ) 1/8 * 1 « 1 <i, 

fiiVH^/H, Vert. Pol. 
fc/«=^ (4-5) 

lV«-*f/*i, Horiz. Pol. 

Here, 2j is the height of the point Xi from the surface, and the function #(£) has been numeri- 
cally calculated by Logan [1959] for the case of a perfect conductor and by Wait and Conda 
[1959] for a wide range of surface impedance. 

Special applications of the supplemental formulas will be as follows: 

Case 1: J m /ij m >l (n>m>2). 
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This is the case where the diffraction angle of every diffracting surface is sufficiently large. 
On using (4.3c) in the formula (3.1), the attenuation coefficient A takes the form 



A={(r n+1 , n +d n + . . . +r 32 +d 2 +r 21 )/kr 1 r 



n+l ,n * n, n-1 



r Z2 r 21 } 1/: 



XMi^a^M^n^a^) . . . M(&, Oa), (4.6a) 

or, when either or both transmitter or receiver is on or near the diffracting surface, the formulas 
(3.2) and (3.3), respectively, give on using (4.5), 



A — Vn+l, n~rd n -\- 



-\-r Z2 +d 2 )lk n l 2 r n +i,nr n ,n-i • • . n 3 r d2 ] 



A={(d n +r nt M _!+ . . . +r S2 + d 2 )/k n l - 3 r n , »_ 



XM(f„, a n ) 



^43^32 j 



M(f 8 , a 3 )flf(&)/(si, a 2 ), (4.6b) 



X/(»»fi, a^tUMf^!, a w _0 . . . M(&, <iz)g(&)f{zi, <h). (4.6c) 

These results just correspond to the so-called multiplication rule in diffraction. 

Another application that may be of interest will be that in which one of the diffracting 
surfaces does not serve as a diffracting obstacle but simply as a reflecting surface, as in figure 24. 
In this case the coupling between the preceding and succeeding surfaces becomes very serious. 
From the leading term of formula (4.3a) we have 

Case 2: £mA?m<— 1 

2_j 1 \?m+L m)t m , lt t m J \£m)t m -L V'm,rn-l) t rn t m _ l 



~T(r m +\, m ) t m+1> o T(r Mj m _i) 0| * TO _! VW(lAv+i, m+ l/r«, m -i) exp [i(U—Q 2 /vl] 
=T(r m +i,m+r m , m -i) t m+h ^_xXexp [i{(U^ 2 -t(r ffl , m .i/(v u +r K(m 4))(k+i/2) 1/3 ^ + i 

-i(^+l t m/(Vl,m+Vm-l))(4l«m-l/2) 1/ ^m-l}] J £m/Vm<> — 1. (4.7) 

When this result is substituted in (3.1), (3.2), and (3.3), they indicate what would happen if 
there were no diffracting surface at all between the m ± 1 th diffracting surfaces: the phase terms of 
B m {r m , m =Fi/(r m+lt m+r m , m _i)}(&ia OT± i/2) 1/3 ^ ±1 are respectively combined with the terms of 
T(£ m± i) t m±l to give the new diffraction angles d' m ± x due to the omission of the mth surface (fig. 
24). The correction to the leading term can be proven to give the reflected wave from the 
mth surface. 

In the same way, when the convergence of the double series 



JLj ■*■ Vm-t-2, m+l) t„ 

tm,t m +l 



l* (£/« + !/ * TO +i^ \ r m-\-l,m)t m ^i,t m J- \&m) t m -*- Vm,m-l)t 



(4.8) 



is poor with respect to both t m and t m +i } we have the leading term (fig. 25) 

J- V m+2,m+l)t m+ 2< ° v m> w-l)o, t m _ 1 { rCl^m+2, m+l^m+1, mTm, m-lf V 'w+2, m+1 ~\^m-\-l, m~T^m,m-V I 



X I { S (f w + i^ /4 , /Wfm+l) +$ (Ue iH \ PJfnd } , (4.9) 



which is valid for ^ w <Cl, r; w+1 ^Cl. Here, 





Figure 24. 77ie situation in the case of eq (4.7). 



Figure 25. The forms of diffracting surfaces for eqs 
(4.8) and (4.9). 
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fim= { U~ (hr m , „,_,) - 1 (k l aJ2y' 3 (k l a m - 1 /2)»H m - 1 }/ v „, 

fm= uPmTPm+l I lvm + 2, ro+l+Pro+1, m) ^m+l,mT^,m-l)/^+l,wVmHw+lT^+l,BJT^,»i-l)] j 
T == I ?w+2, m+lTm, m - l/v'm -41, wi^ot, m-l) Vm+2, b+IT^+1, w) 1 > ©tC. , 



4 fe, 7*)=(-^)V 1+ " 2) * 2 f'tft^-'j fV' 2 ^,. 

YV7T/ J 2 Jn* 2 



(4.10) 



The result of (4.9) corresponds to (4.3a) and can be applied to the formulas (3.1), (3.2), and 
(3.3) for the conditions £ m <CL, £ m+ i<Cl. 

When Zm-i/rjm-i^l and f m+2 /^ m+ 2^>l, the terms of £ m _i and Z m+2 in (4.10) can be neglected 
as in (4.ob), and it follows that 

i m = { (.Vm + l\~Vm) f 'm. m - 1 T^m + l^m + 1, m J {^l^m+2,m-\-ll^Vm-\-l, m\ r m,m-\) V m+2, m+ll^TO+1, /w I Ff», m-l) J > 

^=^V*i/2(l/r w+liW +l/r WtOT _ 1 ), etc, (4.11) 

For the case in which (4.11) can be used, the coupling- between the mth and ra+lth sur- 
faces and the other surfaces is lost, and thus they behave independently in the formulas (3.1), 
(3.2), and (3.3) for the attenuation coefficient A 

5. Diffraction by Two Ridges 

The special case of n=3 and £ 2 ,£ 3 <Cl, for which (4.9) and (4.11) are valid, corresponds to 
the case of two-ridge diffraction, as illustrated in figure 26. A straight application gives the 
result 

Here 



f 2 = { (^2+^3)^21+^32} V^43/2 (7- 32 +r 21 )(r43+r 3 2+r 2 i), 

fs= { (^2+^3) ^3 +^32} V*ir 2 i/2(r 4 3+r32)(r43+r32+/ , 2i), 



ft=^Vfci/2(l/r82 + l/r 21 ),i83=^Vfci/2(l//-32 + l/r43). 



(5.2) 



Equation (5.1) corresponds to the Fresnel integral in the case of a single ridge. Indeed, 
letting r 32 — >0, it follows that 



f 3= fa= (0 8 +0 2 ) -ylkMn/2 (r 43 +r 21 ) , 

and thus, according to the definition of S\z^n) in (4.10) and <%{z) in (4.3a), eq (5.1) becomes 

4=J<f(sV*'*). (5.3) 

Equation (5.1) is given in terms of the function (?(a, b/a) whose precise analytical treatment 
has been studied in another paper [Furutsu, 1955] and also in the appendix of this paper, the 
results of which are briefly described for convenience of computation. It is a many-valued 
function and the following relations are conveniently available for (5.1): 



S\a,b/a)=2e b2 S(a)-cf(a-b/a), a /2 >0, b R <0, 
-<f(-a,6/(-a)), a R <0,b B >0. 
Here a R and b R are the real parts of a and b, respectively. 



(5.4) 



Figure 26. Two ridges and the notations for eqs 
{5.1) and {5.2). 
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Model experiments for two-ridge diffraction have been tried by Decker and his colleagues 
at the National Bureau of Standards. The results of comparison of the theory and the pre- 
liminary experiments are shown in figure 27. Here the theoretical values were obtained by 
the use of a high-speed computer. The agreement is surprisingly good. 
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Figure 27. A few comparisons of experimental and theoretical values for diffraction loss for a double knife- 
path. 

a. Knife edges at variable heights with transmitter and the two edges in line. 

b. Knife edges at equal variable heights, frequency 24.21 Gc/s. 

c. First knife edge fixed, second variable in height, frequency 24.14 Gc/s. 



60 



The result (5.1) was obtained in eq (12.2.15) of the reference [Furutsu, 1956]. Recently 
an equivalent formula was also derived by Millington et al. [1962], and the comparisons with 
the empirical methods of Bullington and Epstein-Peterson were discussed. 



The author expresses his cordial thanks to J. R. Wait for his useful discussions and also 
to M. T. Decker for kindly showing the data of his laboratory experiments on two-ridge 
diffraction. 

Appendix 

From the definition of ${z) in eq (4.3a), 

<?(z) = e z2 -S : (z), J(z)=~e z2 ( Z e- t2 dt. (A.l) 

■yjir Jo 

Here 

n = l 

z 2 2 

Cn= „_1 @n-U Ci=~j= Z. (A. 2) 

lb 2 -\ 7T 

For |s|^>l and |arg(s) |<C37r/4, S; (z) takes the asymptotic forjn 



where 



As)~S g n , M>1, |arg (0)|<3t/4, (A.3) 

// = i 



(A.4) 



n— f 




I 


0*= ^2— 


ff»-i, 


ft =7= 


(O 




V* 



In the same way, from (4. 10), 

S (a, 6/a)=? (tan- 1 £) e« 2+62 -e 62 ?(a) + ?fo, 6/a). (A.5) 

Bere r?(tt) is the same function as in (A.l) and 

— / 2 V o o f " o r {bla)t 2 

£ (a, 6/a)=(4=) ^ + *~ *2«~ ,! ^'^n 

\\7T/ JO Jo 

«f(0, 6/a)= 2 tan- 1 f- (A.6) 

7T 

The function <f(a, 6/a) can be expanded in the absolutely convergent series of 

?(a,W=Si». (A.7) 

Here 

A»=- {(a 2 + 6 2 )A // -i+aJ, a n =- — r b 2 a n - u 

Th Th 2 

A 1=ai =? aft. (A.8) 

For |a 2 +6 2 |*>l and |a/6|<l and |arg(a)|, |arg(6) |<3tt/4, the function *f(a, b/a) 
takes the asymptotic form 
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^{a,bla)~^(^~^iZJ n , |a 2 +& 2 | 1/2 »l, |a/6|<l, |arg (o)|, |arg (ft)|<3r/4. (A.9) 



Here 

' J2n+i = Bj 2n —nAj 2n _ 1 —kn, jo=f?(b), 

32n = Bj2n-l (^' 2j^j2n-2 



{ 



k n =-(n-h)(a 2 +b 2 )- 1 k n . l) h=7r- l/2 b(a 2 +b 2 )-\ 

A=a 2 (a 2 +b 2 )- 2 , B=b 2 {a 2 +b 2 )~\ (A. 10) 

For the case of |a/6|>l, the following relation is applicable: 

£(a, b/a) = S\a) £(fi)-£(fi,a/b). (A.ll) 
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